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Abstract 

The cross section of diffractive deep-inelastic scattering ep — )• eXp is measured, where the 
system X contains at least two jets and the leading final state proton is detected in the HI 
Forward Proton Spectrometer. The measurement is performed for fractional proton longi- 
tudinal momentum loss xjp < 0.1 and covers the range 0.1 < \t\ < 0.7 GeV^ in squared 
four-momentum transfer at the proton vertex and 4 < < 110 GeV^ in photon virtuality. 
The differential cross sections extrapolated to |t| < 1 GeV^ are in agreement with next-to- 
leading order QCD predictions based on diffractive parton distribution functions extracted 
from measurements of inclusive and dijet cross sections in diffractive deep-inelastic scat- 
tering. The data are also compared with leading order Monte Carlo models. 
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1 Introduction 



Diffractive processes such as ep — )■ eXY, where the systems X and Y are separated in rapidity, 
have been studied extensively in deep-inelastic scattering (DIS) at the electroiJll-proton collider 
HERA [Ij-tSJ. Diffractive DIS events can be viewed as resulting from processes in which the 
photon probes a net colour singlet combination of exchanged partons. The photon virtuality 
Q^, the high transverse momentum of jets or a heavy quark mass can provide a hard scale for 
perturbative QCD calculations. For semi-inclusive DIS processes such as ep — j- eXp' the hard 
scattering QCD coUinear factorisation theorem [9] allows the definition of diffractive parton 
distribution functions (DPDFs). The dependence of diffractive DIS on a hard scale can thus 
be treated in a manner similar to the treatment of inclusive DIS, for example through the ap- 
plication of the DGLAP parton evolution equations [fT0] - [T3l . DPDFs have been determined 
from QCD fits to diffractive DIS measurements at HERA [|2l|3l[8l. The inclusive diffractive 
DIS cross section is directly proportional to the sum of the quark DPDFs and constrains the 
gluon DPDF via scaling violations. The production of diffractive hadronic final states contain- 
ing heavy quarks or jets proceeds mainly via boson gluon fusion (BGF) and therefore directly 
constrains the diffractive gluon density [I3][8l. 

In previous analyses at HERA, diffractive DIS events have been selected on the basis of 
the presence of a large rapidity gap (LRG) between system Y, which consists of the outgoing 
proton or its dissociative excitations, and the hadronic final state, system X flSHH. The main 
advantage of the LRG method is its high acceptance for diffractive processes. A complemen- 
tary way to study diffraction is by direct measurement of the outgoing proton, which remains 
intact in elastic interactions. This is achieved by the HI experiment using the Forward Pro- 
ton Spectrometer (EPS) lfT4lfT5l . which is a set of tracking detectors along the proton beam 
line. Despite the low geometrical acceptance of the EPS, this method of selecting diffractive 
events has several advantages. The squared four-momentum transfer at the proton vertex, t, 
can be reconstructed with the EPS, while this is only possible in exclusive final states in the 
LRG case. The EPS method selects events in which the proton scatters elastically, whereas the 
LRG method does not distinguish between the case where the scattered proton remains intact or 
where it dissociates into a system of low mass My. The EPS method also allows measurements 
to be performed at higher values of fractional proton longitudinal momentum loss, xp, than 
possible using the LRG method. 

This paper presents the first measurement of the cross section for the diffractive DIS process 
ep — )■ ejjX'p, with two jets and a leading proton in the final state. The diffractive dijet cross 
sections are compared with next-to-leading order (NLO) QCD predictions based on DPDFs 
from HI [|2l|3l and with leading order (LO) Monte Carlo (MC) simulations based on different 
models. 

The dijet cross sections are measured for two event topologies: for a topology where two jets 
are found in the central pseudorapidity range, labelled as 'two central jets', and for a topology 
where one jet is central and one jet is more forwarc^ labelled as 'one central -i- one forward jet' . 
The universality of DPDFs is studied using events with two central jets. The distributions of 
the proton vertex variables x p and t are compared to those of the inclusive diffractive DIS case. 

'in this paper "electron" is used to denote both electron and positron unless otherwise stated. 
^ The forward direction is defined by the proton beam direction. 
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Figure 1: The leading order boson gluon fusion diagram for dijet production in diffractive DIS. 

This comparison tests the proton vertex factorisation hypothesis which assumes that the DIS 
variable factorise from the four-momentum of the final state proton. The data are also compared 
directly with the LRG measurement of the dijet cross section in diffractive DIS (31 in order to 
test the compatibility of the two experimental techniques. Finally, events with one central and 
one forward jet are used to investigate diffractive DIS in a region of phase space where effects 
beyond DGLAP parton evolution may be enhanced. This topology is not accessible with the 
LRG method since the rapidity gap requirement limits the pseudorapidity of the reconstructed 
jets to the central region. 



2 Kinematics 

Figure [T] illustrates the dominant process for diffractive dijet production in DIS. The incoming 
electron with four-momentum k interacts with the proton with four-momentum P via the ex- 
change of a virtual photon with four-momentum q. The DIS kinematic variables are defined 
as: 

q2= 2^(^_^/)2 x = ^^, y = ^, (1) 
V y V ; , 2p. g' ^ P-k' 

where is the photon virtuality, x is the longitudinal momentum fraction of the proton carried 
by the struck quark and y is the inelasticity of the process. These three variables are related via 
= xys, where s denotes the ep centre-of-mass energy squared. 

The hadronic final state of diffractive events consists of two systems X and Y, separated 
by a gap in rapidity. In general, the system Y is the outgoing proton or one of its low mass 
excitations. In events where the outgoing proton remains intact. My = rup, the mass of the 
proton. The kinematics of diffractive DIS are described by: 

^jp = 5 — ' t = {P-P), p = - — — — — = — , (2) 

q ■ P 2q ■ [P — P') xjp 
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where xp denotes the longitudinal momentum fraction of the proton carried by the colour sin- 
glet exchange, t is the squared four-momentum transfer at the proton vertex and (3 is the frac- 
tional momentum of the diffractive exchange carried by the struck parton. The longitudinal 
momentum fraction of the diffractive exchange carried by the parton entering the hard scatter is 

- q-{P-P'y 
where v is the four-momentum of the parton. 



3 Theoretical Framework and Monte Carlo Models 

Within Regge phenomenology, cross sections at high energies are described by the exchange of 
Regge trajectories. The diffractive cross section is dominated by a trajectory usually called the 
Pomeron {IP). In analyses of HERA data |l2l[3l[8l, diffractive DIS cross sections are interpreted 
assuming 'proton vertex factorisation' which provides a description of diffractive DIS in terms 
of a resolved Pomeron lfT6llT7ll . The QCD factorisation theorem and DGLAP parton evolution 
equations are applied to the dependence of the cross section on and /3, while a Regge inspired 
approach is used to express the dependence on xp and t. 

The resolved Pomeron (RP) model [16] is implemented in the RAPGAP event genera- 
tor [fT8|. RAPGAP implements both a leading Pomeron (iP) trajectory and a sub-leading 
'Reggeon' {FL). In this analysis the DPDF HI 2006 Fit B [0 is used, which employs the 
Owens pion PDFs |fT9l for the partonic content of the Reggeon. The Reggeon contribution is 
significant for xjp > 0.01. Higher order QCD radiation is modelled by parton showers. Pro- 
cesses with a resolved virtual photon are also included, with the photon structure function given 
by the SaS-G 2D LO parameterisation [,20J . 

In the two-gluon Pomeron (TGP) model ll2Tll22l . the diffractive exchange is modelled at 
LO as the interaction of a colourless pair of gluons with a qq or qqg configuration emerging 
from the photon. The model is implemented in the RAPGAP generator. Higher order effects 
are simulated using parton showers. The unintegrated gluon PDF of set AO [|23l is used. 

In the soft colour interaction (SCI) model [|24ll25l . the diffractive exchange is modelled via 
non-diffractive DIS scattering with subsequent colour rearrangement between the partons in the 
final state, which can produce a colour singlet system separated by a large gap in pseudorapidity. 
A refined version of the SCI model which uses a generalised area law (GAL) for the probability 
of having a soft colour interaction [26] is used in this analysis (SCI-i-GAL). Predictions for 
diffractive dijet production within the SCI-i-GAL model are obtained using the leading order 
generator program LEPTO [|27l . Higher order effects are simulated using parton showers [|28l 
|29l . The calculations are based on the CTEQ6L [301 proton PDFs. The probability for a soft 
colour interaction, R, has been tuned to 0.3 to describe the total diffractive dijet cross section 
as measured using the 'two central jets' topology. 

In all three models hadronisation is simulated using the Lund string model fiV\ implemented 
within the PYTHIA program [l321. 
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In this analysis the dijet cross section is also compared to NLO QCD calculations. Assuming 
proton vertex factorisation, NLO QCD predictions for the diffractive partonic dijet cross section 
are calculated in bins of xjp using the NLOJET++ [33 1 program and integrated over the full xjp 
range of the measurement. The renormalisation and factorisation scales are set to /i,. = /i/ = 
^Q^ + (P^)2, where (P^) is the mean of the transverse momenta of the two leading jets in 
the hadronic centre-of-mass frame. In order to estimate the uncertainties of the NLO QCD 
calculations due to missing higher orders, the factorisation scale fij and renormalisation scale 
yUf. are varied simultaneously by factors of 0.5 and 2. The average uncertainty arising from the 
variation of the scale is about 40%. The DPDFs used in the NLO QCD calculations are HI 2006 
Fit B [2J and HI 2007 Jets [3J. The HI 2007 Jets fit is based on the diffractive inclusive and dijet 
data while HI 2006 Fit B is based on inclusive diffractive data only. In order to demonstrate the 
size of the NLO corrections, the QCD calculations are also performed at leading order. 

The NLO QCD partonic cross sections are corrected to the level of stable hadrons by eval- 
uating effects due to initial and final state parton showering, fragmentation and hadronisation. 
The hadronisation corrections are defined in each bin as a ratio of the cross section obtained 
at the level of stable hadrons to the partonic cross sections. Two sets of hadronisation cor- 
rections have been obtained using the RAPGAP generator using two different parton shower 
models: parton showers based on leading logarithm DGLAP splitting functions in leading or- 
der as ifTomSll and parton showers based on the colour dipole model as implemented in ARI- 
ADNE [|34ll . The nominal set of corrections (1 + 5had) is taken as the average of the two sets, 
while the difference between them is considered as the hadronisation uncertainty. The average 
hadronisation corrections are of about 0.9 with an estimated uncertainty of about 7%. Uncer- 
tainties arising due to scale variations and hadronisation corrections are added in quadrature and 
quoted as a total uncertainty on NLO QCD predictions. 

In order to compare with the results of the FPS measurements, NLO and LO QCD predic- 
tions as well as predictions of the RP model are scaled down by a factor of 1.20 [[TSl due to 
the fact that the DPDF sets HI 2006 Fit B and HI 2007 Jets use LRG data which contain a 
proton dissociation contribution. The t-dependence of the iP and ]R fluxes implemented in the 
HI DPDF sets and the RP model are tuned to reproduce the t-dependence measured in inclusive 
diffractive DIS with a leading proton in the final state lfT4l . 



4 Experimental Technique 

The e^p data used in this analysis were collected with the HI detector in the years 2005 to 2007 
and correspond to an integrated luminosity of 156.6 pb^. During this period the HERA collider 
was operated at electron and proton beam energies of E^, = 27.6 GeV and Ep = 920 GeV 
respectively, corresponding to an ep centre-of-mass energy of ^/s = 319 GeV. 

4.1 HI detector 

A detailed description of the HI detector can be found elsewhere [l 35] - [37ll . Here, the components 
most relevant for the presented measurement are described briefly. A right-handed coordinate 
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system is employed with the origin at the nominal interaction point, where the 2;-axis pointing in 
the proton beam or forward direction and the x{y) axis points in the horizontal (vertical) direc- 
tion. The polar angle 9 is measured with respect to the proton beam axis and the pseudorapidity 
is defined as i] = — lntan(6'/2). 

The Central Tracking Detector (CTD), with a polar angle coverage of 20° < 9 < 160°, 
is used to reconstruct the interaction vertex and to measure the momenta of charged parti- 
cles from the curvature of their trajectories in the 1.16 T field provided by a superconducting 
solenoid. Scattered electrons with polar angles in the range 154°<^^g<176° are measured in 
a lead / scintillating-fibre calorimeter, the SpaCal [l37l . The energy resolution is a{E)/E ^ 
1%/ E[GeY] © 1% as determined from the test beam measurement ll38l . A Backward Pro- 
portional Chamber (BPC) in front of the SpaCal is used to measure the electron polar an- 
gle. The finely segmented Liquid Argon (LAr) sampling calorimeter surrounds the track- 
ing system and covers the range in polar angle 4° < 6* < 154° corresponding to a pseudo- 
rapidity range —1.5 < ?7 < 3.4. The LAr calorimeter consists of an electromagnetic section 
with lead as the absorber and a hadronic section with steel as the absorber. The total depth 
varies with 9 between 4.5 and 8 interaction lengths. The energy resolution, determined from 
test beam measurements [ISllllQl, is cr{E)/E ^ ll%/^/E\GeV] © 1% for electrons and 
a{E)/E ^ 50%/^y E[GeV] © 2% for hadrons. The hadronic final state is reconstructed us- 
ing an energy flow algorithm which combines charged particles measured in the CTD with 
information from the SpaCal and LAr calorimeters [|4T|| . 

The luminosity is determined by measuring the rate of the Bethe-Heitler process ep — )■ epj 
detected in a photon detector located at z = —103 m. 

The energy and scattering angle of the leading proton are obtained from track measurements 
in the FPS p2l. Protons scattered at small angles are deflected by the proton beam-line magnets 
into a system of detectors placed within the proton beam pipe inside two movable stations, 
known as Roman Pots. Both Roman Pot stations contain four planes, where each plane consists 
of five layers of scintillating fibres, which together measure two orthogonal coordinates in the 
(x, y) plane. The fibre coordinate planes are sandwiched between planes of scintillator tiles used 
for the trigger. The stations approach the beam horizontally and are positioned at 2; = 61 m and 
z = 80 m. The detectors are sensitive to scattered protons which lose less than 10% of their 
energy in the ep interaction and are scattered through angles below 1 mrad. 

The energy resolution of the FPS is approximately 5 GeV within the measured range. The 
absolute energy scale uncertainty is 1 GeV. The effective resolution in the reconstruction of the 
transverse momentum components of the scattered proton with respect to the incident proton is 
determined to be ~ 50 MeV for and ~ 150 MeV for Py, dominated by the intrinsic transverse 
momentum spread of the proton beam at the interaction point. The scale uncertainties in the 
transverse momentum measurements are 10 MeV for P^ and 30 MeV for Py. Further details 
of the analysis of the FPS resolution and scale uncertainties can be found elsewhere [15J. For 
a leading proton which passes through both FPS stations, the track reconstruction efficiency is 
48% on average. 
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4.2 Kinematic reconstruction 



The inclusive DIS variables Q^, x and the inelasticity y are reconstructed by combining infor- 
mation from the scattered electron and the hadronic final state using the following method [IJ: 

y = ye + yd-yd . Q = , 2(n, io\ , x = — . (4) 

tan^(6';/2) sy 

Here, ye and y^ denote the values of y obtained from the scattered electron only (electron 
method) and from the angles of the electron and the hadronic final state (double angle method), 
respectively [43] . 

The observable xp is reconstructed as: 

xjp = l- E'^/Ep, (5) 

where E'^ is the measured energy of the leading proton in the FPS. The quantity f3 is recon- 
structed as (3 = x/x]p. The squared four- momentum transfer at the proton vertex is recon- 
structed using the transverse momentum Pt of the leading proton measured with the FPS and 
Xip as described above, such that: 

c>2 ^2 ,!->-) 2 

f = f ■ I / . = — (f,) 

1 - Xjp 1- X]p 

where |tmin| is the minimum kinematically accessible value of \t\. The absolute resolution in 
t varies over the measured range from 0.06 GeV^ at |t| = 0.1 GeV^ to 0.17 GeV^ at \t\ = 
0.7GeV2. 

An estimator for the momentum fraction zp is defined at the level of stable hadrons as: 

+ M2, 

zr = , (7) 



xjpys 



where Mjj denotes the invariant mass of the dijet system. The cross sections are studied in 
terms of the DIS variables Q^, (3, zp, the proton vertex variables xjp and t, the jet variables 
and 1], and 

(P^) = ^{Pl, + P^,) , \Av*\ = \vl - r^;\ , |A0*| = 101 - 0^1 , (8) 

where P^^,?]*,^^ and P^2^Vh'p2 transverse momenta, pseudorapidities and azimuthal 
angles of the axes of the leading and next-to-leading jets, respectively, reconstructed in the 
hadronic centre-of-mass frame. The indices 1, 2 stand for the two jets used in the specific 
analyses. 



4.3 Event selection 

The events used in the 'two central jets' and 'one central -i- one forward jet' analyses are trig- 
gered on the basis of a coincidence of a signal in the FPS trigger scintillator tiles and in the 
electromagnetic SpaCal. The trigger efficiency, calculated using events collected with indepen- 
dent triggers, is found to be 99% on average and is independent of kinematic variables. 
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4.3.1 DIS selection 



The selection of DIS events is based on the identification of the scattered electron as the most 
energetic electromagnetic cluster in the SpaCal calorimeter. The energy E'^ and polar angle 
9'^ of the scattered electron are determined from the SpaCal cluster and the interaction vertex 
reconstructed in the CTD. The electron candidate is required to be in range 154° < 9'^ < 176° 
and E'^ > 10 GeV. In order to improve background rejection, an additional requirement on the 
transverse cluster radius, estimated using square root energy weighting [|44ll . of less then 4 cm 
is imposed. 

The reconstructed z coordinate of the event vertex is required to be within ±35 cm of the 
mean position. At least one track originating from the interaction vertex and reconstructed in 
the CTD is required to have a transverse momentum above 0.1 GeV. 

The quantity ~ -^z)' summed over the energies and longitudinal momenta of all re- 

constructed particles including the electron, is required to be between 35 GeV and 70 GeV. For 
neutral current DIS events this quantity is expected to be twice the electron beam energy when 
neglecting detector effects and QED radiation. This requirement is applied to remove radiative 
DIS events and photoproduction background. 

In order to ensure a good detector acceptance the measurement is restricted to the ranges 
4 < g2 < GeV2 and 0.05 <y < 0.7. 

4.3.2 Leading proton selection 

A high FPS acceptance is ensured by requiring the energy of the leading proton E'^ to be greater 
than 90% of the proton beam energy Ep and the horizontal and vertical projections of the trans- 
verse momentum to be in the ranges —0.63 < < —0.27 GeV and \Py\ < 0.8 GeV, respec- 
tively. Additionally, t is restricted to the range 0.1 < \t\ < 0.7 GeV^. 

The quantity ^(.E + P^), summed over all reconstructed particles including the leading pro- 
ton, is required to be below 1880 GeV. For neutral current DIS events this quantity is expected 
to be twice the proton beam energy. This requirement is applied to suppress cases where a DIS 
event reconstructed in the central detector coincides with background in the FPS, for example 
due to interactions between off- momentum protons from the beam halo with residual gas within 
the beampipe. 

Previous diffractive dijet DIS measurements |l3l|4l[6]| and DPDF fits [|2l[3l[8l have been per- 
formed for |tmin| < N < 1 GeV^. To compare with these results, the cross sections are 
extrapolated to the range |tmin| < 1^1 < 1 GeV^ using the t dependence measured in inclusive 
diffractive DIS with a leading proton in the final state lfT4l . 

4.3.3 Jet selection 

Reconstructed hadronic final state objects are used as input to the longitudinally invariant fcy jet 
algorithm [|45l using thepr recombination scheme with ajet radius of 1.0 as implemented in the 
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FastJet package [|46ll . The jet finding algorithm is applied in the photon-proton centre-of-mass 
system (7*]? frame). The jet variables in the 7*p frame are denoted by a asterisk. 

In the 'two central jets' analysis, the requirements are P^j^ > 5 GeV and > 4 GeV 
for the leading and next-to-leading jet, respectively. Asymmetric cuts are placed on the jet 
transverse momenta to restrict the phase space to a region where NLO calculations are reliable. 
The axes of the jets are required to lie within the pseudorapidity range — 1 < r^i 2 < 2.5 in 
the laboratory frame. The selected event topology is similar to that in the LRG dijet data used 
in the DPDF fits [|3l[8]|. This data selection is used for testing the proton vertex factorisation 
hypothesis and the DPDFs in processes with a leading proton in the final state. 



Selection 


two central jets 


one central + one forward jet 


DIS 


4 < g2 < iioGeV^ 




0.05 < y < 0.7 


Leading Proton 




xjp < 0.1 






t\ < 1 GeV2 




Pf 1 > 5 GeV 


Pt,c, Ptj > 3-5 GeV 


Jets 


P^ 2 > 4 GeV 


Mjj > 12 GeV 




-1 < r/1,2 < 2.5 


-1 <r]c< 2.5 






l<Vf < 2.8, r^j > r/c 



Table 1: Phase space of the diffractive dijet FPS measurements. 

The selection of the 'one central + one forward jet' topology is motivated by the study of 
diffractive DIS processes in a phase space where deviations from DGLAP parton evolution 
may be present. The requirement of a forward jet suppresses the parton pt ordering which is 
assumed by DGLAP evolution. At least one central jet with —1 < rjc < 2.5 and one forward jet 
with 1 < rif < 2.8, where rff > r]c, are required with P^ > 3.5 GeV. In addition, the invariant 
mass of the central-forward jet system is required to be larger than 12 GeV to avoid the phase 
space region in which NLO QCD calculations are unreliable. 

The selection criteria for the two analyses are summarised in table [H The 'two central jets' 
data sample contains 581 events and the 'one central -1- one forward jet' data sample contains 
309 events. 



5 Corrections to the Data and Cross Section Determination 

5.1 Background subtraction 

The selected data samples contain background events arising from random coincidences of non- 
diffractive DIS events, with off-momentum beam-halo protons producing a signal in the FPS. 



11 



The beam-halo background contribution is estimated statistically by combining the quantity 
'^{E + Pz) summed over all reconstructed particles in the central detector in DIS events (with- 
out the requirement of a track in the FPS) with the quantity XI (-^ + ^z) beam-halo protons 
from randomly triggered events. The J2iE + Pz) spectra for leading proton and beam-halo 
DIS events for both dijet event topologies are shown in figure |2l The background distribution 
is normalised to the FPS DIS data distribution in the range '^{E + P^) > 1880 GeV where 
the beam-halo background dominates. The ratio of signal to background depends on the signal 
cross section and is found to be considerably larger than in the inclusive diffractive DIS pro- 
cesses measured with the FPS detector [15J. After the selection cut + ^z) < 1880 GeV 
the remaining background amounts on average to about 5%. The background is determined and 
subtracted bin-by-bin using this method. 

5.2 Detector simulation 

Monte Carlo simulations are used to correct the data for the effects of detector acceptance, 
inefficiencies, migrations between measurement intervals due to finite resolution and QED ra- 
diation. The response of the HI detector is simulated in detail using the GEANT3 program ll47l 
and the events are passed through the same analysis chain as is used for the data. The reaction 
ep —7- eXp is simulated with the RAPGAP program [18] using the RP model and the DPDF set 
HI 2006 Fit B as described in section [3l QED radiative effects are simulated using the HERA- 
CLES [|48l program within the RAPGAP event generator. In the 'two central jets' analysis the 
r]2 distribution of the Monte Carlo simulation is reweighted in order to describe the experimen- 
tal data. A similar procedure is applied to the r]j distribution in the 'one central -i- one forward 
jet' sample. More details of the analysis can be found elsewhere [ 49J . 

A comparison of the FPS data and the RAPGAP simulation is presented in figure [3] for the 
variables xp and \t \ reconstructed with the FPS detector. The contributions of light quarks (uds) 
to ]P and ]R exchanges and of charm quarks to ]P exchange are also shown in the logiQ(xjp) 
distribution. Figure |4] presents the data and the Monte Carlo distributions of the variables 
P^i, ^iid Z]p for the 'two central jets' sample and of the variables {P^),T]f and zp for 
the 'one central -i- one forward jet' topology. For this comparison zjp is reconstructed from the 
scattered electron and the hadronic final state in the HI detector. The MC simulation reproduces 
the data within the experimental systematic uncertainties. 

5.3 Cross section determination 

In order to account for migration and smearing effects and to evaluates the dijet cross sections 
at the level of stable hadrons, matrix unfolding of the reconstructed data is performed [50] . The 
resolution and acceptance of the HI detector is reflected in the unfolding matrix A which relates 
reconstructed variables y^ec with variables on the level of stable hadrons Xtrue via the formula 
Axtrue = Vrec- The matrix A, obtained for each measured distribution using the RAPGAP 
simulation, is constructed within an enlarged phase space in order to take into account possible 
migrations from outside of the measured kinematic range. The following sources of migrations 
to the analysis phase space are considered: migrations from low Q'^, from low y, from large xp. 
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from low Pt jets, from the single jet topology, fulfilling the Pt requirements for the leading jet 
as given table [H and in case of the 'one central + one forward jet' analysis from large r]f. In 
order to treat the contamination of the measurement by these migrations correctly the analysis 
is performed in an extended phase space which includes side-bins in y^ec and xtme for each of 
the migration sources listed above. 

The unfolded true distribution on the level of stable hadrons is obtained from the measured 
one by minimising a function defined as 

X^ = XI + r\l = l/2(y;cc - AXtrnefY~\yrec ' Aftrue) + r\l (9) 

where Xa ^ measure of a deviation of Axtrue from the data bins ^rec- The matrix V is the 
covariance matrix of the data, based on the statistical uncertainties. In order to avoid statistical 
fluctuations, the regularisation term Xl implemented into the function and defined as 
Xl = (a^truc)^- The regularisation parameter r is tuned in order to minimise the bin-to-bin 
correlations of the covariance matrix V. Further details of the unfolding method can be found 

iniEDISIl 

The Bom level cross section is calculated in each bin i according to the formula: 

(Tiiep ^ ejjX'p) = j{l + 5rad) (10) 

where Xi is the number of background subtracted events as obtained with the unfolding proce- 
dure described above, C is the total integrated luminosity and (1 + Srad) are the QED radiative 
corrections which amount to about 5% on average. The differential cross sections are obtained 
by dividing by the bin width. 



6 Systematic Uncertainties on the Measured Cross Sections 

The systematic uncertainties are implemented into the response matrix A and propagated through 
the unfolding procedure. They are considered from the following sources listed below. 

• The uncertainties on the leading proton energy and on the horizontal and vertical projec- 
tions of the proton transverse momentum are 1 GeV, 10 MeV and 30 MeV, respectively 
(section l4~T1) . The corresponding average uncertainties on the cross section measurements 
are 0.5%, 5.3% and 2.2%. The dominant uncertainty originates from the FPS acceptance 
variation as a function of the leading proton transverse momentum in the horizontal pro- 
jection. The above uncertainties result from the run-by-run variations of the incoming 
proton beam angle and of the FPS detector positions relative to the proton beam, as well 
as from the imperfect knowledge of the HERA beam magnet optics. 

• The uncertainties of the measurements of the scattered electron energy E'^ (1%) and angle 
6'g (1 mrad) on the SpaCal calorimeter lead to an average systematic uncertainty of the 
cross section of 1.5% and 2.8%, respectively. 
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• The systematic uncertainty arising from the hadronic final state reconstruction is deter- 
mined by varying the energy scale of the hadronic final state by ±2% as obtained using a 
dedicated calibration [i53il . The 2% uncertainty of the calibration is confirmed by studies 
in the region of low jet transverse momenta and low photon virtuality. This source leads 
to an average uncertainty of the cross section measurements of 6.2% for production of 
two central jets and 9.5% for production of one central and one forward jet. 

• The model dependence of the acceptance and migration corrections is estimated by vary- 
ing the shapes of the distributions in the kinematic variables (P^), r]2, rjj, xp, (3 and 

in the RAPGAP simulation within the constraints imposed on those distributions by the 
presented data. The rjl and rj*^ reweightings are varied within the errors of the parame- 
ters of the reweighting function, which amount up to a factor 4. The (P^) distribution is 
reweighted by (Pf the xp distribution by (l/xip)=^° °^, the /3 distribution by 
and (1 — and the Q"^ distribution by log((5^)^°'^. For the 'two central jets' se- 

lection the largest uncertainty is introduced by the rfl reweighting (4%), followed by (3 
(2.7%), while the reweights in xjp, (P^) and result in an overall uncertainty of 2.3%. 
The uncertainties for the 'one central -i- one forward jet' topology are 12.8% for the r^^ 
reweighting, followed by (P^) (2.1%), while the reweights in xp, (3 and Q"^ result in an 
overall uncertainty of 1.8%. 

• Reweighting the t distribution by e^* results in a normalisation uncertainty of 4.2% for 
the extrapolation in t from the measured range of 0.1 < \t\ < 0.7 GeV^ to the region 
l^minl < 1^1 < 1 GeV^ covered by the LRG data [3]. The uncertainty arising from the 
t reweighting within the FPS acceptance range of 0.1 < \t\ < 0.7 GeV^ is on average 
1.4%. 

The following uncertainties are considered to influence the normalisation of all measured 
cross sections in a correlated way: 

• Two sources of systematics related to the background subtraction are taken into account: 
the energy scale uncertainty and the limited statistics in the data sample without the 
^{E-\-p;^) cut. Firstly, the beam-halo spectrum is shifted within the quoted uncertainties 
of the hadronic energy scale and proton energy scale. Secondly, the normalisation of the 
background spectrum is shifted by 1 ± 1/ ^/N^, where Ai'bkg is the number of events in 
the FPS data sample in the range + -^z) > 1880 GeV. The uncertainties from these 
two sources are combined in quadrature. The uncertainty of the proton beam-halo back- 
ground is considered as a normalisation error and found to be 3.5% for the production of 
two central jets and 1.5% for the production of one central and one forward jet. 

• A normalisation uncertainty of 1% is attributed to the trigger efficiencies, evaluated using 
event samples obtained with independent triggers. 

• The uncertainty in the FPS track reconstruction efficiency results in a normalisation un- 
certainty of 2%. 

• A normalisation uncertainty of 3.7% arises from the luminosity measurement. 
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The systematic errors shown in the figures are obtained by adding in quadrature all the contri- 
butions except for the normalisation uncertainties, leading to an average uncertainty of 11% for 
'two central jets' and 17% for 'one central + one forward jet'. The overall normalisation un- 
certainty of the cross section measurement obtained by adding in quadrature all normalisation 
uncertainties is 7% for 'two central jets' and 6.2% for 'one central -i- one forward jet'. The cross 
section measurement in t has a normalisation uncertainty 4.6%. 

7 Results 

The ep cross section for diffractive production of two central jets and one central -i- one forward 
jet, integrated over the full measured kinematic range (table[T]), is given in table [21 together with 
the predictions obtained with NLO QCD calculations. 





two central jets 
a[ph] 


one central -i- one forward jet 
^[pb] 


Data 


254 ±20 (stat.)±27 (syst.) 


150 ± 19 (stat.) ± 26 (syst.) 


NLO QCD 
HI 2006 FitB 
HI 2007 Jets 


270 tlf (scale) ± 16 (hadr.) 
257 til (scale) ± 22 (hadr.) 


148 tlf (scale) ± 6 (hadr.) 
128 tlr (scale) ± 7 (hadr.) 



Table 2: Total cross section for the 'two central jets' and 'one central -i- one forward jet' samples 
compared to the NLO QCD calculations. 

Within the uncertainties, both cross sections are well described by the NLO QCD calcula- 
tions. 

The measured differential cross sections are presented in tables [BUS] and figures [5lfT4l The 
tables also include the full covariance matrices of the experimental uncertainties. The quoted 
differential cross sections are averaged over the intervals specified in the tables [BHSl 

7.1 Differential cross section for the production of two central jets 

The measured differential cross sections are shown as a function of Q^, y, logiQ^xjp) and zp 
in figure [51 The calculations obtained with the DPDF sets HI 2006 Fit B and HI 2007 Jets are 
presented as well as the ratio R of the calculations to the data. Within the uncertainties, the 
normalisation and shape of the cross sections are reasonably well described by the NLO QCD 
predictions. Since dijet production is directly sensitive to the gluon DPDF, the measured cross 
sections confirm the normalisation and shape of the gluon DPDFs extracted from the NLO QCD 
fits to diffractive inclusive and dijet cross sections measured using the LRG method [[2l[3l. 

In figure [6l the differential cross sections are shown as a function of ^ and | Ar^*]. Within 
the errors, NLO QCD predictions describe the data. A slight deviation of the theory from the 



15 



data is observed for jets with a small separation in pseudorapidity |Ar7*|. For the differential 
cross section in P-^i, the LO QCD contribution is calculated as well using the DPDF set HI 
2007 Jets and is observed to underestimate the measured cross section by a factor of about 2. 

Figure |7] shows a comparison of the differential cross sections in Q'^, y, log^o(^jp) ^^'^ 
with MC models based on the leading-logarithm approximation and parton showers. The ratios 
of the measured cross sections to the MC predictions show that the RP model gives a good 
description of the shape, but underestimates the dijet cross section by a factor of 1.5. For this 
comparison the reweighting with respect to the rfl distribution specified in section 15.21 is not 
applied to the RP model. Since the IP and ]R fluxes which determine the xjp dependence in the 
RP model has been tuned to the inclusive diffractive DIS LRG data [2J the good agreement in 
shape of the RP model with the dijet data supports the hypothesis of the proton vertex factori- 
sation. Both the SCI-i-GAL and TGP models fail to describe the data. The SCI-i-GAL model 
predicts harder spectra in and zp and a softer spectrum in log^^^xp) than are seen in the 
data. It should be noted that the probability of soft colour interactions and hence the normal- 
isation of diffractive processes in the SCI-i-GAL model is adjusted to the measured dijet cross 
section. The TGP model is in agreement with the data only at low xjp but underestimates the 
data significantly at larger x jp where the missing sub-leading contributions are expected to be 
large. 

Figure [8] shows the differential cross sections in P^^ and |Ar7*| for the data and the MC 
models. The shapes of these distributions are again well described by the RP model. Although 
the SCI-I-GAL model is not able to describe the differential cross sections as a function of the 
diffractive kinematic variables xjp and zjp and of the DIS kinematic variable this model 
reproduces reasonably well the measurements as a function of the jet variables P^ and | Ar7*|. 

None of the LO Monte Carlo models are able to describe all features of the measured dif- 
ferential cross sections. The best shape description in all cases is provided by the RP model. 
However, this model is a factor of 1.5 below the data in normalisation. The TGP and SCI-i-GAL 
models fail to describe the shape of the differential cross sections. 

The differential cross section in \t\ shown in figure [9^ is fit using an exponential form 
exp(i?t) motivated by Regge phenomenology. An iterative procedure is used to determine 
the slope parameter B, where bin centre corrections are applied to the differential cross sec- 
tion in t using the value of B extracted from the previous fit iteration. The final fit results in 
B = 5.89 ± 0.50 (exp.) GeV'^, where the experimental uncertainty is defined as the quadratic 
sum of the statistical and systematic uncertainties and the full covariance matrix is taken into ac- 
count in the fit. As shown in figure[9j5, this t-slope parameter is consistent within the errors with 
the t-slope measured in inclusive diffractive DIS with a leading proton in the final state IfTSl at 
the same value of xjp. The consistency of the measured t dependence with that for the inclusive 
diffractive DIS cross sections supports the validity of the proton vertex factorisation hypothesis. 

The cross section for the production of two central jets can be compared with the diffractive 
dijet measurement obtained using the LRG technique |[3l. The LRG measurement includes 
proton dissociation to states Y with masses My < 1.6 GeV. To correct for the contributions 
of proton dissociation processes, the LRG dijet data are scaled down by a factor of 1.20, taken 
from the diffractive inclusive DIS measurement [[T5l . To compare to the results of the LRG 
method, dijet events are selected in the same kinematic range. The DIS and jet variables Q^, y. 
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P-f 1 and r]i^2 are restricted to the ranges 4 < < 80 GeV^ 0.1 < y < 0.7, ^ > 5.5 GeV, 
and —1 < r^i 2 < 2, respectively. The results are presented in figureflOl The comparison shows 
consistency of the results within the experimental errors. Compared to the LRG measurement, 
the phase space of the present analysis extends to xp values that are a factor of three larger. 

7.2 Differential cross section for the production of one central + one for- 
ward jet 

Figure [TT] shows the differential cross sections for the production of 'one central + one forward 
jet' as a function of | Ar]* \,r]f and the mean transverse momentum of the forward and central jets 
(P^) together with the expectations from the NLO QCD. Within the errors, the measured data 
are described by NLO QCD predictions. In order to test the predictions in a wider kinematic 
range, the ?7j distribution of the forward jet shown in figure [TT] is extended down to a minimum 
value of —0.6 where the prediction overshoots the data. LO QCD calculations, performed using 
the DPDF set HI 2007 Jets underestimate the measured cross section by a factor of about 2.5. 

The differential cross sections measured as a function of zp, log^Q(/3) and |A0*| are pre- 
sented in figure [121 The data are well described by the NLO QCD predictions. In the BFKL 
approach [54-56], additional gluons can be emitted in the gap between the two jets, leading to a 
de-correlation in azimuthal angle | A</)*|. The observed agreement between the measured cross 
sections and NLO DGLAP predictions in this distribution shows no evidence for such an effect 
in the kinematic region accessible in this analysis. 

Figure [T3] presents the differential cross sections for the production of 'one central -i- one 
forward jet' as a function of the variables (P^), \A.r]*\ and rjj. The RP model is a factor of 2.2 
below the data which is a larger discrepancy in normalisation than that observed in the 'two 
central jets' sample. A similar trend is seen for the LO QCD contributions in the two samples. 
The normalisation of the SCI-i-GAL model, tuned to 'two central jets', agrees with the cross 
section for 'one central -i- one forward jet' . The shapes of the distributions are reasonably well 
described by both the RP and SCI-i-GAL models. 

The differential cross sections in zjp, \ogiQ{/3) and | A0*| are shown in figure[l4] The shapes 
of all distributions are well described only by the RP model. As for the case of the 'two central 
jets' the SCI-I-GAL model is not able to describe the distributions of the diffractive kinematic 
variables but it well reproducing the shape of the | A0* | distribution. The TGP model completely 
fails again to describe the zp spectrum. 



8 Summary 

Integrated and differential cross sections are measured for dijet production in the diffractive 
DIS process ep ejjX'p. In the process studied, the scattered proton carries at least 90% of 
the incoming proton momentum and is measured in the HI Forward Proton Spectrometer. The 
presented results are compatible with the previous measurements based on the LRG method and 
explore a new domain at large xp. 
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Dijet cross sections are measured for an event topology with two jets produced in the central 
pseudorapidity region, where DGLAP parton evolution mechanism is expected to dominates, 
and for a topology with one jet in the central region and one jet in the forward region, where 
effects of non-DGLAP parton evolution may be observed. NLO QCD predictions based on the 
DGLAP approach and using DPDFs extracted from inclusive diffraction measurements describe 
the dijet cross sections within the errors for both event topologies, supporting the universality 
of DPDFs. The measured t-slope of the dijet cross section is consistent within uncertainties 
with the value measured in inclusive diffractive DIS with a leading proton in the final state. 
This confirms the validity of the proton vertex factorisation hypothesis for dijet production in 
diffractive DIS. 

The measured cross sections are compared with predictions from Monte Carlo models based 
on leading order matrix elements and parton showers. The Resolved Pomeron model describes 
the shape of the cross sections well, but is too low in normalisation. This suggests that contri- 
butions from higher order processes are expected to be sizable in this approach. The SCI+GAL 
model is able to reproduce the normalisation of the cross section for both dijet topologies pre- 
sented after tuning the model to the 'two central jets' data. The dependence of the diffractive 
dijet cross section on xjp and zp is able to distinguish between the models. The SCI-i-GAL 
and Two Gluon Pomeron models fail to describe the shape of the distributions of the diffractive 
variables, while the Resolved Pomeron model describes the shape of these distributions well. 
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Table 3: Bin averaged hadron level differential cross sections for the production of two central jets in diffractive DIS as a function of , y, 
log^o(^J') ^iid zp. The total {5tot), statistical {5stat) and systematic {dgys) uncertainties and the correlation coefficients p of the cross section 
covariance matrix defined in section |53] are given together with the changes of the cross sections due to a +lcr variation of the various 
systematic error sources described in section [6l the electromagnetic energy scale {5eie), the scattering angle of the electron {5e), the leading 
proton energy Ep (6ep), the proton transverse momentum components (SpJ and Py (Spy), the reweighting of the simulation in ri2 (5^2 ) 
and xjp (Sxjp) the hadronic energy scale (5Ehad), the reweighting of the simulation in (3 (5^), (5q2) and P^ (Sp*). All uncertainties are 
given in per cent. The normalisation uncertainty of 7% is not included. The hadronisation correction factors (1 + Shad) applied to the NLO 
calculations and the associated uncertainty are shown in the last column. 
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Table 4: Bin averaged hadron level differential cross sections for the production of two central jets in diffractive DIS as a function of Pt^i, 
\Ar]* I and \t\. The normalisation uncertainties of 4.6% for the differential cross section in \t\ and 7% for other cross sections are not included. 
For details see table [3l 
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Table 5: Bin averaged hadron level differential cross sections for the production of one central and one forward jet in diffractive DIS as a 
function of {P^), I Vf^ ^ip^ logio(/^) 1^0*1- The normalisation uncertainty of 6.2% is not included. For more details see tabled 
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Figure 2: The distribution of ^(i? + P^) for FPS DIS events (points) and for beam-halo DIS 
events (shaded histogram). 
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Figure 3: The distributions of the variables xp (a) and |t| (b) reconstructed using the FPS 
(points) for events with two central jets. The beam-halo background is subtracted from the data. 
The RAPGAP Monte Carlo simulation, reweighted to describe the distribution, is shown as 
a histogram. Contributions from sub-processes are illustrated in the xjp distribution as areas 
filled with different colours. 
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Figure 4: The distributions of the variables P^i, |A?7*| and zjp for events with two central jets 
and of the variables {Pt), Vf for events with one central and one forward jet (histogram 

with the error bars). The beam-halo background is subtracted from the data. 
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Figure 5: The differential cross section for the production of two central jets shown as a function 
of Q^, y, logiQ{x]p) and zp. The inner error bars represent the statistical errors. The outer error 
bars indicate the statistical and systematic errors added in quadrature. NLO QCD predictions 
based on the DPDF set HI 2007 Jets, corrected to the level of stable hadrons, are shown as a 
solid line and a dark green band indicating the hadronisation uncertainties and light green band 
indicating the hadronisation and scale uncertainties added in quadrature. LO QCD predictions 
based on the same DPDF set are shown as a dotted line. The NLO calculations based on the 
DPDF set HI 2006 Fit B with applied hadronisation corrections are shown as a thick line. R 
denotes the ratio of the measured cross sections and QCD predictions to the nominal values of 
the measured cross sections. The total normalisation error of 7.0% is not shown. 
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Figure 6: The differential cross section for production of two central jets shown as a function 
of I and \Ari* \ . For more details see figure [51 
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Figure 7: The differential cross section for the production of two central jets shown as a function 
of Q^, y, \ogiQ{x]p) and zjp. The inner error bars represent the statistical errors. The outer error 
bars indicate the statistical and systematic errors added in quadrature. The RP and SCI+GAL 
models are shown as solid and dotted lines, respectively. R denotes the ratio of the measured 
cross sections and MC model predictions to the nominal values of the measured cross sections. 
The total normalisation error of 7.0% is not shown. 
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Figure 8: The differential cross section for production of two central jets shown as a function 
of ^ and I A?7*|. For more details see figure |71 
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Figure 9: The differential cross section for production of two central jets shown as a function of 
if: (a), the corresponding t- slope (circle) shown as a function of xp (b). The result is compared 
to the HI inclusive diffractive DIS data (triangles) [|T5l . The error bars indicate the statistical 
and systematic errors added in quadrature. 
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Figure 10: The differential cross section for the production of two central jets in the phase space 
of the LRG measurement [S]] as described the text in section 17.11 The cross section is shown 
as a function of log^Q(a;jp). The inner error bars represent the statistical errors. The outer error 
bars indicate the statistical and systematic errors added in quadrature. The published LRG dijet 
data are scaled down by a factor of 1.20 to correct for the proton dissociation contribution are 
shown as open circles with the error bars indicating the statistical and systematic errors added 
in quadrature. 
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Figure 1 1 : The differential cross section for the production of one central and one forward 
jet shown as a function of the mean transverse momentum of two jets (-Rf)' Vf- 
The inner error bars represent the statistical errors. The outer error bars indicate the statistical 
and systematic errors added in quadrature. NLO QCD predictions based on the DPDF set HI 
2007 Jets, corrected to the level of stable hadrons, are shown as a line with a dark green band 
indicating the hadronisation error and light green band indicating the hadronisation and scale 
errors added in quadrature. LO predictions based on the same DPDF set are shown as a dotted 
line. The NLO calculations based on the DPDF set HI 2006 Fit B with applied hadronisation 
corrections is shown as a thick line. R denotes the ratio of the measured cross sections and QCD 
predictions to the nominal values of the measured cross sections. The total normalisation error 
of 6.2% is not shown. 
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Figure 12: The differential cross section for production of one central and one forward jet shown 
as a function of zp, logio(/3) and |A0*|. For more details see figure [TTl 
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Figure 13: The differential cross section for production of one central and one forward jet 
shown as a function of the mean transverse momentum of two jets (P^), |A?7*| and r]f. The 
inner error bars represent the statistical errors. The outer error bars indicate the statistical and 
systematic errors added in quadrature. The RP and the SCI+GAL models are shown as solid 
and dotted lines, respectively. R denotes the ratio of the measured cross sections and MC model 
predictions to the nominal values of the measured cross sections. The total normalisation error 
of 6.2% is not shown. 
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Figure 14: The differential cross section for production of one central and one forward jet shown 
as a function of zp, logio(/3) and |A0*|. For more details see figure [T3l 



35 



